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Abstract In this paper, we study some fourth order singular critical equa- 
tions of Lichnerowicz type involving the Paneitz-Branson operator, and we 
prove existence and non existence results under given assumptions. 

1 Introduction 

During the last years there have been effective studies of conformal operators 
and their relative invariants due to their application in geometry or mathe- 
matical physics. For instance the Yamabe problem played an essential role in 
the evolution of the analytical and geometrical tools also it was with crucial 
importance for the study of the Einstein-Hilbert functional without forget- 
ting the input in relativity for the study of the conformal Einstein constraint 
equations (see [9], [8]). And In 1983, Paneitz [20J introduced a conformally 
fourth order operator defined on 4-dimensional Riemannian manifolds Bran- 
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son [3j generalized the definition to n— dimensional Riemannian manifolds, 
n > 5. He introduced another geometric quantity that defines another con- 
formal invariant for n > 3 that is the Q — curvature that behaves in a very 
similar way to the scalar curvature. And its variation after a conformal 
change involves a fourth order operator. That is the Paneitz-Branson Op- 
erator. One can think about the Q — curvature and the Branson-Paneitz 
operator like the scalar curvature and the conformal laplacian. There was 
a lot of published work concerning prescribing the Q— curvature where one 
can notice that the conditions that we get are similar to the scalar curva- 
ture one modulo some technical assumptions (see |13|.|14j.[2]). One of the 
issues that we meet while dealing with this operator, is the fact that there 
is no maximum principle, thus getting good and effective estimate is not as 
easy as for the conformal Laplacian. Many authors have studied the posi- 
tivity and coercivity of the Paneitz, one can consult [23] or [16] for example. 
As interaction with mathematical physics we can see the work of Choquet- 
Bruhat in [9] with the Conformal Laplacian, where they study the scalar 
field equation that leads to a Lichnerowicz type semi-linear PDE. In this 
work we attempt to study another action functional as a proposal for a rela- 
tivestic model since it is conformally invariant and we will see a scalar-field 
perturbation of it. The study of such functional leads to the resolution of 
a Lichnerowicz type equation but it is a fourth order one, with the Paneitz 
operator as a differential part. So in this paper we will investigate the ex- 
istence of positive solutions under some assumption that we will mention 
later to that equation and also we will give a non-existence result. 
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2 Preliminaries and Motivations 



Let (M,g) be a n— dimensional closed compact manifold with n > 3, recall 
that if R g is the scalar curvature then under the conformal change g = u n ~ 2 g, 
one gets the following relation relation the new curvature with the old one : 



(n — 2) n+2 

gU + 4(n-l) gU = gU "~ 2 ' 



Let —L g u = — A g u + 4^~l^Rg u this operator is called the conformal Lapla- 
cian to see more of its property one could check [18] . Similarly if we consider 
the following quantity which is the Branson Q — curvature introduced in 
[3], defined by 



Q : 
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Then after a conformal change g = u n ~ 4 g of the metric, one gets 



n+4 



Q U n-4 = P U 



(2) 



where 



(n-2) 2 + 4 _ 4 _. \ , \ n-4^ 
i?5 -Ric \du\-\ —Qu. 



Pu := A 2 r ,u — div . . , ... 

9 I \2(n-2)(n-l) n-2 



We will set Po its differential part, that is 

Po-ii = - div \\ ^ n ~ 2 } ' +4 —Pic | du | , 

5 U 2(n-2) («- 1) y n-2 i /' 

One can see that if 5 is an Einstein metric then P is with constant coefficient, 



Pu = Aiu + 



n 1 -2n- A 
2n(n - 1) 



PA n + 



(n - 4) (n 2 - 4) 
16n(n — l) 2 



P^n 



and satisfies the maximum principle since it can be written as a product of 
two second order operator satisfying the maximum principle. 

Remark that the natural space to work on for the prescribed scalar cur- 
vature problem is H 1 (M) and thus ([!]) is a critical semi-linear problem. Also 
the natural space for prescribing the Branson Q-curvature is H 2 (M) and 
again ([2]) is a critical problem since we are in the borderline of the Sobolev 
embeddings. 

There was an extensive work concerning <Q to find a metric with con- 
stant scalar curvature, which is a critical point of the Einstein-Hilbert func- 
tional 



restricted to the conformal Class of a given metric. Same thing can be 
applied to the functional 



For instance one could check that the functional is Reimannian and 
Einstein metrics are critical points of this functional. 

There have been many proposal in relativity to replace the Hilbert- 
Einstein total curvature functional with a conformally invariant functional, 
like for instance the case of Bach relativity where the functional is replaced 
by f M \Cg\ 2 dv g and C is the Weyl tensor (see [3]). In this case we will 
consider another proposal consisting of the total Paneitz-Branson curvature 
Fq. Therefore one can think about a scalar field perturbation of the previous 
one, that is 
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this functional was studied in [10] for the case Hilbert-Einstein functional 
under conformal change where the authors try to solve a conformal con- 
straint for the Einstein scalar field equation also in [15] where the author 
studies the problem from a variational point of view. 

Now if we take a closer look to this functional, one can see that if we 
restrict it to the conformal class of g one have 
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Fib ( u n ~ 4 g ) = — / uPu — a n \V „ty\ 2 u 2 du q . 
V / a n J M 

Where a n = rk ^, therefore, the associated Euler-Lagrange equation to this 
problem is 

Pg^u = Pu - a n \V g ^\u = [Q - |Vg*| )«»-*, 

that is 

. 2 \ ( 7\ I--. T |2\ "+ 4 



P u + a n [Q- \V g % ) u= [Q- |Vg*| J u«-4 
and the constant 

«>0,neC oo (M) a n / „ 2n \ — 

{jM Un - 4 ) 

is a conformal invariant. 

Let us recall the following result about the coercivity of the Paneitz- 
Branson operator. 

Theorem 2.1 ([M])- Let (M,g) be a closed Riemannian manifold of di- 
mension at least 6. If the Yamabe invariant of g is non-negative, then with 
respect to any conformal metric of positive scalar curvature Po has a non- 
negative first eigenvalue and kerP $ = {constant} . The last statement also 
holds in dimension five, provided we assume the Yamabe metric has positive 
Q- curvature. 



Proposition 2.2. Under the assumptions of the previous theorem, then 
the sign of V [g, is the sign of the first eigenvalue of the Operator Pg^ 
for every g in the same conformal class. 

Proof. Assume V [g, VP] > 0, then if we take (pi the first eigenvalue as a test 
function, one gets 

Ai|M! = E(pi)>7>[0,tt]||pi|& 

thus Ai > 0. Now if we assume that V [g, \T>] = then using the same 
argument we get that Ai > , but if Ai > then using Sobolev inequalities 
we get V [g, VP] > which is not the case, thus Ai = 0. 

Now if V [g, \t] < then there exist a function v, such that E{v) < 
thus we get Ai < 0. □ 

From now on we will assume that the Yamabe and the Paneitz invariants 
are positive and P is positive, therefore we guaranty that the operator P g ^ 
is coercive and satisfies the maximum principle and up to a conformal change 
we can assume that = Q — |Vt/>|^ id positive on M. 

And if we follow the procedure of the Authors in |10| to find the Einstein- 
scalar field conformal constraint equation one gets a a Lichnerowicz type 
problem but of fourth order of the following form : 

P 9 ,*u = $L - B{x)«*-i 
u > 

where 2" = ^rj, A and B two smooth functions. Therefore, the object of 
the rest of this paper is to investigate the existence of positive solutions to 
problem of the following form : 

[ P g ^u=^-B(x)u<i 
u > 
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where p > 1 and 1 < q < 2" — 1. 

3 Existence Via heat flow 

Let £ be a Banach space with norm || — \\ E . E is partially ordered by a 

o 

closed cone P C X, and we assume that it has non-empty interior P. We 
define also P = P — {0}. The element of P are called positive and element 
of — P are called negative. 

Now, if u, v € E we will use the following notations to distinguish how 
they are comparable : 

u <v \iv — u £ P 
< u<v\iv — u£P ■ 

o 

u <^Lv it v — u £ P 

v 

A map /, we set D(f) C E its domain. Now a map / : D(f) — > E is said 
order preserving, if for every u,v G D(f) such that u < v then f(u) < /(u). 
And we say that 

lim fix) = +oo 

x — >+oo 

if for every u £ P there exist x £ P such that /(f) > u, for every w > x. 
And finally we define the set 

[u, v ] = {w £ E; u < w < v} 

and sometimes if needed for a set D C E, 

[u, v] D = {w £ D; u < w < v} . 

An OBS is said normal if there exist 5 > such that for every u < v in 

E, 

\\u\\ E < S \\v\\ E . 
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Theorem 3.1 (Krein-Rutman) . Let E be a total ordered Banach space and 
T : E — > E a compact order preserving linear operator, then r(T) is an 
eigenvalue with eigenvector u £ P and if in adition we assume that T is 
strongly order preserving (That is Tu >> if u > 0) then r(T) > and is 
a simple eigenvalue with positive eingen vector. 

Let us consider the following problem 

4iu + Au = F{u) 
< dt ' (3) 

u(0) = n 

o 

where F : P — > X is a C 1 map , A is a densely defined compact resolvent 

o 

positive operator and uq G Pfl D(A). 
Theorem 3.2. Assume that 

lim F(x) = — oo 

x — >+oo 

and 

lim Fix) = +oo 

and for every bounded set K there exist a a constant A such that F + XI is 
order preserving in K, then the problem admits a positive solution. 



In the applications we can know more about the solution and we will 
deal with that further in this paper. 

Proof. First remark that there exist u\ and U2 such that 

^ui+Aui < F(m) 
ui(0) < u 
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and 

^u 2 + Au 2 > F(u 2 ) 

o 

In fact u\ and u 2 could be chosen of the form se where e G P and s > 0. 
set 

K = {u\ < u < u 2 } 

then if is a bounded set, so there exist A > so that F + XI is order 
preserving on K so let A and F denote respectively A + XI and F + XI. 

Now let us construct the following sequence : u 1 being the unique solu- 
tion of 

+ Au = F(ui) 
and u k+1 is the unique solution of 

^u + Au = F(u k ) 

By induction one can easily show that the sequence (n fc ) is monotone 
non-decreasing and u k G K, \/k > 1. Let us show the first step, that is 
u 1 > U\. 

First using the assumptions on A we have the existence of a compact 
positive semi-group S(t), generated by A (see |21j). So we have 

ui = S(t)ui(0) + J S{t - s) (J^Ui + Am 

< S(t)u + [ S{t- s)F( Ul (s))ds 

Jo 

< u l 

Now since A has compact resolvent and K is bounded we can extract for 
fixed time a subsequence that we will call also (u fc ) such that S(t)u k (s) 



converges to S(t)u(s), thus by writing 

u k+1 = S(t)u + f S(t-s)F(u k ), 
one can see that u satisfies 



o 



u = S(t)u + [ S(t-s)F(u 
Jo 



And this gives a positive solution to 

Now notice that u(t) is bounded in D(A) thus there exists a sequence 
(tfe)t, going to infinity such that u(t).) converges to some u, and in fact 
the convergence occurs in D(A). Thus knowing that lim f Q S(t — s)xds = 
{—A)~ 1 x we get by passing to the limit that u is a solution of the steady- 
state problem. □ 

Now we will consider a problem of the form 

U t + Pg^U = f(x,u) 

u(0) = u 

where P is the Paneitz-Branson operator and / : M x — > R is C 1 such 

that lim x >+OQ f(x,t) = — oo uniformly on x and lim^. ^Qf(x,t) = +00 

uniformly on x, Then for every uq smooth and positive, there exist a positive 
solution to Q moreover there exist a sequence tk going to infinity such that 
u(tk) — > u a solution of the steady-state problem. 

One also can write the problem as an integral equation using the posi- 
tivity preserving flow like in 



Corollary 3.3. Take A and B two positive functions defined on M and 
consider the singular problem : 

Pu=^- B(x)ui 
u > 



10 



then using Theorem (3.2) we have the existence of a solution more than that, 
it is the unique solution. 

Remark that in this case we can take q > 2$ — 1 since we do not need 
the compact or continuous embedding in L p spaces. 

Now, if we suppose that B is just non-negative, we can show indeed that 
even in that case we have a solution. 

Corollary 3.4. Take A > and B > 0, two smooth functions defined on 
M and consider the singular problem : 

[ Pg,^U = - B(x)li q 

u>0 

where q < 2$ — 1, then it has a unique solution. 

Let u £ be the solution obtained by Corollary (3.3), of 

( ' Pg^u = ^ - B £ {x)u* 
u>0 

V 

where B £ = B + e. First remark that u £ is uniformly bounded from 
bellow (it is by construction of the sub and super solution in the proof of 
Theorem (3.2)). 

So 

/ uP 9 ^u =[*&-[ {B(x)+e)u^ 

JM J M U £ J M 

< [ *f = c 

Jm o 

where 5 is the uniform lower bound of u £ . Therefore (u £ ) £ is bounded in 
H 2 (M) and if q + 1 < 2", we have u £ — > u in L? and weakly in H 2 (M) 
and L 2 \M). 
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So if we take 92 G C°°(M), we have a weak solution which we can show 
using the regularity theory that is is indeed smooth. 

<pP g>i> u £ = [ ^-cp - [ (B{x) + e) u%<p 

M JM Ue JM 

so by letting e — > we get that 

f A(x) f 
ipP g ^u= / — j-y- / B(x)u q ip 

M JM UF JM 

so u is a weak solution and using elliptic regularity we get the fact that it 
is indeed a smooth one. 

For the uniqueness, if we consider two smooth positive solutions u and 
v of then w = u — v satisfies : 

= — C(x)(u — v) = — C{x)w 

where C(x) is a non-negative function that we get from the mean value 
theorem, therefore using the maximum principle we get the desired result. 
As an improvement of the previous result we have : 

Theorem 3.5. Assume that B + is non-zero then problem |5p has at least 
one positive solution if the following inequality is satisfied 

max U^rVi^ P+q < ^ FT, ( 6 ) 



M 



£+1 , . 3=1 ' 
q— 1 ^ P+1 1 / p+1 ^ P+1 

p+l) + U-l 



where Ai and (pi are the first eigenvalue and eigenfunction of Pgrf, respec- 
tively. 

first let u be a solution of 

P g ,ipu = ^-B+(x) U o 
u>0 
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In fact since we are going to use this process another time let us give the 
picture and the idea behind : 

Consider a convex function positive / : K + — > R and, so for it to in- 
tersect a line L passing through the origin its slope should be greater than 
the one of the unique tangent to the graph of passing through the origin as 
shown in the following figure : 




and the condition to find the slope at zero is by solving 

/(*) 



t 

So if we take (pi the first eigenfunction of P, we get 



fit)- (7) 



t p^-^ + B{x)t^\ = a m -^| + i?(z)tv? 

> tXnp! - - B-{x)t q <p\. 



And here we can see that in fact we are comparing t — > t\\ipi and t 



, Vi Jp + B (x)t q ip1 which is convex, thus using the previous remark we can see 



A{x) 

that the inequality ^ insures that we are in the same situation as ?? and 
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thus there exist to > such that to<pi is a super-solution to ([5]) therefore, 
using the classical monotone iteration method we get a positive solution. 



3.1 Further investigations and existence results 

Here we investigate the case where B < 0. The coercivity assumption implies 
that 



defines a norm equivalent to the H 2 (M) norm. So we will use that norm 
instead of the usual one. Also we take the best Sobolev constant with 
respect to this norm, that is is the est constant satisfying 



Remark that for B < this condition still work, but let us try to find an- 
other condition that works in a weaker setting. We will rewrite the problem 
as 



and B here is taken to be positive (in fact one get a similar result if B has a 
negative part up to a small modification to the assumption in the following 
theorem) . 

For the regularity issues we refer to [5] and [5], there one can find the 
necessary regularity and bootstrapping argument to deal with it. 

Theorem 3.6. Assume that P is strongly positive (that is it satisfies the 
strong maximum principle). If there exist a function ip > in H 2 (M) such 
that 





Si> < \\u\\tp 





(9) 
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then problem |§p has at least one positive smooth solution. 
In fact we will compute an exact value of C, That is 

c = s ; a W a ( («-')&■-!) ) (10) 

Befor Starting the proof let us state the following lemma which appears to 
be helpful in our situation. 

Lemma 3.7. Let E, E\, E2 be three C 1 functional on a Banach space X. 
Assume that E\(0) = and limE(t<p) = —00. and E2 > 0. Then If E\ has 
the montainpass geometry around zero, (that is there exist r > such that 
5 = 'udggr r \ E\ > 0) and there exist u £ B(0,r) such that E2(u) < 5, the 
functional E has a Palais-Smale sequence. 

Proof of Lemma. Here is is easy to see that if we concider the set 

T = {7 : [0, 1] — > X such that 7(0) = u and 7(1) = tp} 

then we get a Palais-Smale sequence at the level 

c= inf max E(n{ [0,1])) 

Since each curve crosses dB(0,r), then c > max(E(u), E(tip)), and thus we 
have a mountain pass geometry. □ 

Proof. In fact let ipbea positive function such that ||y|L = 1 and the energy 
functional 

E (u) = - \\u\\l ■, + — I r-—[ Bui +l . 

Clearly the functional E\ defined by 
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M 



-1 g+1 

has the mountain pass geometry and in fact if ro = H-BHJ^ 1 S, 2(<J_1) then 



inf Ei_(u) = II-BII?! 1 " ~ ' X 



8+1 / Q — 1 

- - We S ^ ' 

ue9B(0,r ) V 2 



And therefore the inequality ^ is exactly saying that that top satisfies the 
assumption of Lemma (3.7) for to < ro and thus we have the existence of 
to < < t 2 such that, 

max(E(t (p),E(t 2 (p)) < E(r <p) 

And in fact we can apply the lemma for the following approximated energy 
functional 

for e > and small. Remark that we have uniform convergence of t — > 
E E (tcp) to t — > E(t(p), on every compact of R. Therefore there exist £o > 
such that for every < £ < £q one have 

max(E £o (to(p),E £Q (t 2 (p)) < max.(E £ (t <p) , E £ (t 2 ip)) < E £ (r (p) < E(r tp). 

(11) 

Therefore if we take 

r = {7 : [0, 1] — > H 2 (M) such that 7 (0) = t ip; 7 (1) = t 2 ip} 

we have a Palais-smale sequence at the level 

c £ = inf max E £ (u) 
7er«e7([o,i]) 



notice that from (11) 



inf Ei{u) < c £ < E(r tp) 

u£dB(0,r ) 
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and thus c £ is uniformly bounded. Let us call that Palais-Smale sequence 
(u e k )k, it satisfies then 

E £ (u%) — > c £ and E' £ (u%) — > in H~ 2 as k — > oo. 

Thus The following holds 



0(\\u%\ 



(9-1) 



/ £ ll 2 



q+l 
p-1 



A 



Therefor 



0(1), 



which implies the boundedness of in H 2 (M) and thus the existwwence 
of u £ G F 2 (M) such that 

u| — 1 u £ weakly in H 2 
u k — ^ n £ strongly in L 2 , 
ii| — u e weakly in L 2 

so take rj € C°°(M), the previous assertion gives that 

/ vPg^Ue = [ AUeV p -i + / Bujri 

JM JM ( . / + \2\^ JM 

[e+{ut) ) 
thus u £ is a weak solution to the problem 



Au+ 



Pg,ip u e — 

(e + (ut) 2 ) 

hence u £ is smooth and positive. 

First, assume that (jpj — > 0, then 



JW+Buj 



(14) 



(q + 1) E £ {u £ ) - (E' £ {u £ ), u £ ) = {q+ 1)) c £ 
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therefore from (12) we get 



A 



M 



e + (u E 



\2\ 2 



r < Ci 



(15) 



< Co 



where C%, C2 are constants independant of e. Thus we can extract a subse- 
quence of (u e ) £ that we will call (n £fc ) £fe so that 

tt weakly in H 2 (M) 



strongly in L S (M) for 1 < s < 2$ 



u a.e on M. 



Thus using Fatou's lemma in (j 1 5 p we get 

1 



/ 

J M 



(16) 



Assume now that there exist Xk — > x such that u Ek (xk) — > 0. Then using 
the integral representation we get 



u £k (x k )> G(x k ,y)B (y) u q £k {y)dy 

JM 

where G is the Green's function of the operator P g1 b. Taking k — > we get 

G(x,y)B(y) U i(y)dy = 0, 

M 



thus u = which is impossible because of (16), therfor u £ is uniformly 
bounded from below. 



So now we can pass to the weak limit in (14) to get. 



M 



■qPg^u= f —+ f Bu q rj, for every r] £ C°°(M), 

JM uP JM 



hence, since u is positively bounded from below, we get a smooth positive 
solution to 

p g^ u = ~£ + Buq - 
vP 
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If p — 1 = q + 1 (and that is the case of the Lichnerowicz Equation), to find a 

uniform bound on J M — ^rr, we use the fact that ||u e |L is uniformly 

( £ +(« £ ) 2 )^ 

bounded, and Sobolev embedding to get a uniform bound on J M Bu% +1 and 
thus we get the desired bound. □ 

Corollary 3.8. Under the assumption of the previous theorem, we have the 
existence of another positive solution 

Proof. If we take a look at the inequality ^ we notice that it is open, 
that is iw we pertube B a small perturbation, we still get the same ex- 
istence result. So let us call ub the solution corresponding to B. Then 
using a comparison principle, we get ub- £ < ub+s and they are a pair of 
sub and super-solution to the problem Q, therefore we have the existence 
of a solution u to the problem, and to guaranty that ub 7^ u we use a 
degree theory argument since every positive smooth solution is in the set 
A = {itG C 4 ' a (M) ; ^ < -u < C} for C > large enough and uniform. □ 

Corollary 3.9. There exist a constant C = C(n,M,Q^) > such that, if 
Pgib is strongly positive and 

3n-4 f 

(maxB)^ / A < C 

JM 

the Paneitz- Lichnerowicz Problem admits at least one positive solution. 

For the proof of this corollary, we just take tp = 1 in ^ , and the Sobolev 
embedding 

H 2 (M) «-»■ L 2> (M) . 
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4 Non existence Result 



Theorem 4.1. Assume that A, B > 0, then if 



(/ 



(p+q)(q-3) 

g v \ q(p+q-2) 
Ap+i B i+p 



Q 

p+1 



(p+q)2 \ 

q(p+q-2) \ 



+ 



p+1 

p+q-2 



> 



J 



M 



then the problem does not posses any positive smooth solution, where 

Q-|Wf 



Proof. Let u be a positive solution, then the follwong holds : 

A 



+ / Bu q 



/ Q^u 

JM 



using the fact that 



a; 



9-1 

i \ q 



Also 



/ 



q p 
Ap+i Bi+p < 



A \ p+q 



Bu q 



M 



q+p 



therefore if we set X = f Bu q , one gets 



X + 



(/ 

\Jm 



Ap+i B i+p 



p+q 



Which is equivalent to say that 



X 

Therefore if 



\Jm 



p+q 
q p \ q 



Ap+i B i+p 



q-i 
q 



a 



(p+q)(q-3) ( 

g p \ q(p+q-2) 
Ap+i Bi+p 

\l J 



(p+q)2 

q(p+q-2) 



+ 



P+1 

q-1 \ p+q-2 
p+1 



> 



(/ 

\Jm 



then there is no smooth positive solution to the problem. 



□ 
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5 Conclusion 



As a conclusion of the previous existence and non-existence result, we can 
set for the sake of simplicity, A = 1, B = A £ 1, and we get the following 
corollary if we consider the following problem 

PgrfU =— + (17) 

Corollary 5.1. If P g ^ is strongly positive, then there exist a constant X* > 
such that 

i) Problem has no positive smooth solution if X > A*. 

ii) Problem has at least one positive solution if X < A*. 
Moreover we have the following estimate 

q-1 

2»-„-l r \-H)\f+l (e+aVa-ll / r, _ 1 \ vvla-2 9<.P+9-2) 

pq+q-2 
1 

q-1 



2^-g-i / n -4 r A (P P+1 (p+q)(q-i) ( a - 1 \ pp+i-2 

Vol(M) Z-C\—-JqA j < A < Vol(M) p!+q=2~~ ^ i__-J HQ^ii 



where C is the constant (10). 
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